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, We present new results on the cosmological evolution of the near-infrared galaxy 

luminosity function, derived from the analysis of a new sample of ~ 22, 000 Kab< 22.5 
galaxies selected over an area of 0.6 square degrees from the Early Data Release of the 
UKIDSS Ultra Deep Survey (UDS). Our study has exploited the multi- wavelength 
coverage of the UDS field provided by the new UKIDSS WFCAM K and J-band 
M \ imaging, the Subaru/XMM-Newton Deep Survey and the Spitzer-SWIKE Survey. The 

unique combination of large area and depth provided by this new survey minimises 
the complicating effect of cosmic variance and has allowed us, for the first time, to 
trace the evolution of the brightest sources out to z ~ 2 with good statistical accuracy. 

In agreement with previous studies we find that the characteristic luminosity of 
the near-infrared luminosity function brightens by ~ 1 magnitude between z — and 
z ~ 2, while the total density decreases by a factor ~ 2. Using the rest-frame (U — B) 
colour to split the sample into red and blue galaxies, we confirm the classic luminosity- 
dependent colour bimodality at z < 1. However, the strength of the colour bimodality 
is found to be a decreasing function of redshift, and seems to disappear by z > 1.5. 
Due to the large size of our sample we are able to investigate the differing cosmological 
evolution of the red and blue galaxy populations. It is found that the space density 
of the brightest red galaxies (Mk < —23) stays approximately constant with redshift, 
and that these sources dominate the bright-end of the luminosity function at redshifts 
z < 1. In contrast, the brightening of the characteristic luminosity and mild decrease 
in space density displayed by the blue galaxy population leads them to dominate the 
bright-end of the luminosity function at redshifts z > 1. 
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1 INTRODUCTION 
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The potential of deep near-infrared observations to measure 
the cosmological evolution of the stellar mass of galaxies has 
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long been understood (e.g. Lilly & Longair 1984; Dunlop et 
al. 1989; Glazebrook et al. 1995; Cowie et al. 1996). Now, 
with the advent of wide-field near-infrared detectors, this 
potential has begun to be realized, and several deep infrared- 
based studies supported by spectroscopic and/or photomet- 
ric redshift information have recently been completed (e.g. 
Cimatti et al. 2002; Pozzetti et al. 2003; Dickinson et al. 
2003; Drory et al. 2003; Fontana et al. 2004; Glazebrook et 
al. 2004; Caputi et al. 2006; Drory et al. 2005). 

A key aim of this work is the determination of the cos- 
mological evolution of the if-band galaxy luminosity func- 
tion (LF), the present-day form of which is now reasonably 
well established (Cole et al. 2001; Kochanek et al. 2001). 
The most recent studies have been deep enough, and have 
been armed with sufficient redshift information to attempt 
to trace this evolution out to z = 2 — 3. The results indicate 
a substantial brightening in characteristic magnitude from 
z = 0uptoz = 2- 3 along with a simultaneous decrease in 
galaxy number density (see e.g. Caputi et al. 2006; Saracco 
et al. 2006). A similar combination of luminosity and density 
evolution has also been observed in the optical LF (e.g. Lilly 
et al 1995; Poli et al. 2003; Gabasch et al. 2004; Giallongo 
et al. 2005). 

However, even these most recent deep near-infrared sur- 
veys have still only covered relatively small areas of sky 
(from a few arcmin 2 in the Hubble Ultra Deep Field to 
several hundred arcmin 2 in the Great Observatories Ori- 
gins Deep Survey (GOODS)) and are hence existing re- 
sults on the evolving LF are still potentially heavily af- 
fected by cosmic variance. In this paper we exploit the 
order-of-magnitude improvement in areal coverage (0.6 sq. 
deg.) provided by the the UKIRT Infrared Deep Sky Survey 
(UKIDSS; Lawrence et al. 2006) Ultra Deep Survey (UDS) 
to constrain the evolution of the near-infrared LF and the 
colour bimodality up to z > 2. 

The rest of this paper is structured as follows. In Sec- 
tion 2 we summarize the main properties of the datasets 
used in this work. Then, in Section 3 we describe in detail 
how our new K-band galaxy sample has been constructed 
and cleaned, and describe the method we have developed 
for the estimation of redshifts from the available multi-band 
photometry. In Section 4 we assess and discuss the reliabil- 
ity of our photometric redshift estimates and in Section 5 
we derive the near-infrared LF. The evolution of the colour 
bimodality is presented in Section 6 while in Section 7 we de- 
rive separate LFs for the populations of red and blue galax- 
ies. Our conclusions are presented in Section 8. 

Throughout this paper all magnitudes are given in the 
AB system (Oke & Gunn 1983) and we have adopted a 
background cosmology with Qm = 0.3, Qa = 0.7 and 
H = 70 km s _1 . 



2 THE DATA-SETS 

The data exploited in this paper cover a large range in wave- 
lengths from the optical to the mid-infrared and have been 
obtained by combining three different surveys. The optical 
data were taken with the Subaru Telescope as part of the 
Subaru/XMM-Newton Deep Survey (Sekiguchi et al. 2005), 
the near-infrared data form part of the early-data release 
(Dye et al. 2006) of the UKIRT Infrared Deep Sky Survey 



(Lawrence et al. 2006) and mid-infrared data have been ob- 
tained with the Spitzer satellite as part of the Spitzer wide- 
area infrared extragalactic (SWIRE) survey (Lonsdale et al. 
2003; 2004). In this Section the main properties of these 
three surveys are briefly outlined. 

2.1 The Subaru/XMM-Newton Deep Survey 

The Subaru/XMM-Newton Deep Survey (SXDS) is a multi- 
wavelength survey covering an area of ~ 1.3 square degrees. 
This area, centred on RA=02:18:00 and Dec=-05:00:00 
(J2000), benefits from deep optical imaging undertaken with 
Suprime-Cam (Miyazaki et al. 2002) on Subaru. The 5 
over-lapping Suprime-Cam pointings provide broad band 
photometry in the BVRi'z' filters to typical 5<r depths of 
B = 27.5, V = 26.7, R = 27.0, i' = 26.8 and z = 25.9 
(within a 2-arcsec diameter aperture). 

The X-ray observations of this field were obtained 
with the XMM-Newton satellite and comprise 400 ksec 
exposures in 7 contiguous fields to a depth of ~ 
few 10~ 15 erg cm -2 s _1 . The field also has deep 1.4 GHz ra- 
dio observations from the VLA (Simpson et al. 2006a) and 
850pim sub-mm observations from SCUBA as part of the 
SHADES survey (Mortier et al. 2005; Coppin et al. 2006). 

2.2 The UKIDSS Ultra Deep Survey 

The central region of the SXDS is being observed at near- 
infrared wavelengths with the Wide-Field Camera (WF- 
CAM - Casali et al. 2006, in preparation) on UKIRT, as 
part of the UKIDSS Ultra Deep Survey (UDS). The UDS is 
the deepest of the five surveys that constitute the UKIRT 
Infrared Deep Sky Survey (UKIDSS - Lawrence et al. 2006) 
and ultimately aims to cover 0.8 square degrees to a depth 
Jab = 26, Hab = 25.4 and Kab = 25 over a seven year 
period. AB magnitudes in the UDS have been obtained by 
using calibrations from Hewett et al. (2006). For the present 
study the first twelve hours of J, K imaging, which con- 
stitute the early data release (EDR - Dye et al. 2006) of 
the UDS, have been used. These data reach 5a limits of 
Jab = 22.5 and Kab = 22.5 within a 2-arcsec diameter 
aperture. The details of the stacking and catalogue extrac- 
tion are fully described in Foucaud et al. (2007). 

2.3 The SWIRE survey 

The Spitzer wide-area infrared extragalactic survey (SWIRE 
- Lonsdale et al. 2003; 2004) is a large Spitzer Legacy pro- 
gram which covers seven high-latitude fields with a total 
area of ~ 65 sq. deg. In this paper we exploit the Spitzer ob- 
servations covering 9.1 sq. deg. in the XMM-LSS field which 
cover the SXDS field. The source catalogue based on the 
data release 2 (Surace et al. 2005) includes sources with de- 
tections at IRAC 3.6/im and 4.5/im wavelengths to a depth 
of ~ lO^Jy SNR ~ 10 in both bands ( ~ 21.4 in AB mag- 
nitudes) . 

3 THE UDS-GALAXY SAMPLE 

Exploiting the multi-wavelength coverage and the large area 
provided by the combination of the surveys described in the 
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Figure 1. Photometric redshift plotted versus spectroscopic redshift for ~ 1100 sources in the UDS sample with available spectroscopic 
redshift determinations (Yamada ct al. 2005; Akiyama et al. in preparation; Smail et al. in preparation; Simpson et al. 2006a). Left panel 
shows the comparison for sources with Kab < 22 while the comparison for fainter sources with Kab > 22 is presented in the right 
panel. The small insets in both panels show the distribution of the Az/(l + z). The agreement is remarkably good with a mean value of 
Az/(1 + z) of 0.006, and a standard deviation a = 0.04, excluding the clear outliers (see text for details). 



previous Section, we aim to derive a complete sample of K- 
band selected galaxies, with photometrically estimated red- 
shifts. Only sources with Kab< 22.5 have been included in 
the sample, which is the 5a magnitude limit of the EDR data 
(see Section 2.2). The completeness at the limiting magni- 
tude has been estimated to be above 70% with a fraction 
of spurious detections less than 1% (Foucaud et al. 2007). 
In order to have an homogeneous multi- wavelength coverage 
for each object we use optical and near-infrared photometry 
computed in a 2" diameter apertures PSF corrected to to- 
tal. The consistent image quality of optical and near-infrared 
data (0.7 < FWHM < 0.8 arcsec) means that the differen- 
tial aperture corrections between bands are small (< 0.1 
mag). The Spitzer SWIRE photometry is in 2.8" diameter 
apertures PSF corrected to total (see Surace et al. 2005). 

Due to a small shift of the UDS field centre during the 
first year of observations, the current J, K photometry in 
the EDR is uniform only over an area of 0.6 square degrees. 
Therefore, in order to obtain homogeneous coverage over 
all the optical and near-IR spectral range, we have limited 
our sample to include only sources in the central 0.6 square 
degree area of the UDS field. 

Many of cross-talk in the J and K-band images have 
been masked in the process of image stacking and removed 
from the catalogue (see Foucaud et al. 2007). Moreover, 
given the surface density of optical sources and the 2" aper- 
ture diameters used to compute the photometry we expect 
each cross talk to have < 0.03 counterparts in the optical 
catalogues. All the sources with only J and K but not optical 
counterpart have been excluded from the final catalogue. 

Excluding the saturated and masked regions (due to 
bright stars) our final sample comprises 25,564 sources with 
Kab< 22.5. 



All of the sources in this sample have optical and near- 
infrared photometry, and a cross-match with the SWIRE 
catalogue also provided 3.6pim and 4.5/xm fluxes for 14,879 
sources (i.e. more than half of the sample). 



3.1 Photometric redshifts 

The photometric redshift for each galaxy in our UDS/SXDS 
sample has been computed by fitting the observed photom- 
etry (BVRi'z'JK and 3.6/um and 4.5/um when available) 
with both synthetic and empirical galaxy templates. The 
fitting procedure is undertaken with a code which is largely 
based on the public package HYPERZ (Bolzonella, Miralles 
& Pello 2000). The stellar population synthesis models of 
Bruzual & Chariot (2003) were used to produce the syn- 
thetic templates, assuming a Salpeter initial mass function 
(IMF) with a lower and upper mass cutoff of 0.1 and 100 
Mq respectively. We used a variety of star formation his- 
tories; instantaneous burst and exponentially declining star 
formation with e-folding times 0.3 < r(Gyr) < 15, all with 
a fixed solar metallicity. The dust reddening has been taken 
into account by following the obscuration law of Calzetti et 
al. (2000) within the range < A v < 2. We also added 
a prescription for the Lyman series absorption due to the 
HI clouds in the inter galactic medium, according to Madau 
(1995). Finally, at each redshift only models with age less 
than the age of the Universe at that redshift have been al- 
lowed. 

As empirical templates we used the four Coleman, Wu 
and Weedman (CWW) observed spectra: Ellipticals, Sbc, 
Scd and Irr (Coleman, Wu & Weedman 1980). In addition 
we used three average templates obtained by the K20 sur- 
vey (Cimatti et al. 2002); red passive early-type galaxies, 
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intermediate galaxies with emission lines but red continuum 
indices and blue emission- line galaxies (Mignoli et al. 2005). 
An observed starburst SED from Kinney et al. (1996) was 
also included. All of the empirical templates were extended 
into the ultraviolet and infrared wavelengths by using the 
most appropriate synthetic model from Bruzual & Chariot 
(2003). 

Finally, a quasar template was also included to allow 
the photometric redshift code to identify objects with type- 
I AGN features. For this purpose we used the compos- 
ite quasar spectra from the FIRST Bright Quasar Survey 
(Brotherton et al. 2001). 

3.2 Sources of contamination 

Since our primary goal is to construct a sample of Tf-band 
selected galaxies, it is important to clean our catalogue of 
stars and AGNs. To isolate stars we first used the Sextrac- 
tor (Bertin & Arnout 1996) stellarity parameter as mea- 
sured in the i' band in the SXDS images. However, due to 
the ~ 0.8" FWHM resolution of the available ground-based 
imaging ~ 6500 objects were consistent with being unre- 
solved point sources. Therefore, we combined the stellarity 
information with the position of the stars in the BzK colour- 
colour diagram (Daddi et al. 2004). As pointed out by Daddi 
et al., stars have colours that are clearly separated from the 
region occupied by galaxies and can be efficiently isolated 
with the criterion (z' — K) < 0.3(B — z') — 0.5. By applying 
this criterion, coupled with the requirement of stellarity pa- 
rameter > 0.8, we isolated 1990 stars (see Lane et al. 2007, 
for a full discussion of the features displayed by the BzK di- 
agram for our new UDS/SXDS sample). It is worth noting 
that the vast majority of these sources show an unacceptable 
value of the x 2 when fitted with the galaxy templates. 

By using the quasar template in the SED fitting pro- 
cess we also isolated 909 sources for which the photometry 
is dominated by AGN light and therefore better fitted by 
the AGN template. These objects are characterised by a 
blue UV-optical continuum which is virtually featureless in 
the broad-band photometry. Therefore, for these sources the 
photometric redshift technique (which relies on the identi- 
fication of continuum features, such as the 4000A break in 
the case of galaxies) is very hard to apply. The uncertain- 
ties related to the adopted quasar template from the FIRST 
survey (Brotherton et al. 2001) could in principle introduce 
some bias in the fitting procedure, specially due to the con- 
tribution of galaxy light to the composite template (see e.g. 
Glikman et al. 2006 and Maddox & Hewett 2006 for a de- 
tailed discussion) . However, this is not a major issue for the 
present work since the quasar template has only been used 
to isolate very blue featureless AGN therefore excluded from 
the final galaxy sample. The uncertainties related to the 
quasar template will however only affect < 1% of the total 
galaxy sample. Moreover the surface density of such objects 
in the UDS/SXDS field (~ 0.4/arcmin 2 ) is in agreement 
with that obtained in the GOODS/CDFS. In fact, Caputi 
et al. (2006), identified nearly 60 AGN in their .^-selected 
sample down to Kab< 22.5 over an area of 131 arcmin 2 . 

Finally we excluded from the final sample 624 sources 
for which the SED fitting procedure provided a bad \ 2 ■ 
Specifically, given the number of data points and degrees 
of freedom, we reject all fits for which the value of the \ 2 
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Figure 2. A plot of observed K-band magnitude versus photo- 
metric redshift for the galaxies in the new UDS EDR sample. 
The dashed line corresponds to the empirical K — z relation for 
radio galaxies as obtained by Willott et al. (2003), which approx- 
imately coincides with the passive evolution of a 3L* starburst 
formed at redshift zj = 10. The solid line corresponds to the 
passive evolution of a 10L* starburst. 

falls outside the 99.7% confidence interval of the correspond- 
ing x 2 distribution. A visual inspection of the SED fits for 
these sources revealed that their photometry is contami- 
nated, mainly due to an erroneous cross-match between the 
optical and NIR catalogues or source blending. 

3.3 The final UDS-galaxy sample 

In summary, the final sample of galaxies in the UDS/SXDS 
field (hereafter the UDS-galaxy sample) consists of 22, 041 
sources with Kab< 22.5, over an area of 0.6 square degrees, 
excluding stars, AGN and sources with contaminated pho- 
tometry. All the sources in the final UDS-galaxy sample have 
optical + near-IR photometry and for nearly half of them 
we also possess 3.6/im and 4.5/xm detections. 



4 THE REDSHIFT DISTRIBUTION OF THE 
UDS-GALAXY SAMPLE 

The SED fitting procedure provides two vital pieces of infor- 
mation for each galaxy: the best fitting galaxy template and 
a photometric redshift. In this Section we check the reliabil- 
ity of our photometric redshift estimates and compare the 
redshift distribution of sources in the UDS-galaxy sample 
with that derived by previous studies in different fields. 

4.1 Spectroscopic redshifts 

Figure 1 shows the comparison of our photometric redshift 
estimates with the spectroscopic values for ~ 1100 galax- 
ies in the UDS field spanning < z < 6 (Yamada et al. 
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2005; Akiyama et al. in preparation; Smail et al. in prepa- 
ration; Simpson et al. 2006a). We observe a good agree- 
ment between estimated and real redshifts in the vast ma- 
jority of the cases with only < 3% of clear outliers. The 
A«/(l + z) = (zspoct — 2 phot )/(l + 2 S pcct) has a mean of 
—0.003 with a standard deviation a = 0.1. By excluding 
the clear outliers (with |Az/(l + z)\ > 0.2) wc obtain a 
mean of 0.006 with a a reduced to 0.04. This accuracy is 
comparable to the best available from other surveys such 
as GOODS and COSMOS (Caputi et al. 2006; Grazian et 
al. 2006; Mobasher et al. 2007). The residual small focus- 
ing of the photo-z technique on the scale of Sz < 0.05 does 
not affect any of our results since our binning in redshift to 
compute LF and colour bimodality is at least of 0.25. 

The accuracy of our photometric redshift estimates is 
preserved also at faint K-band magnitudes. In fact, ~ 10% of 
the sources with available spectra have magnitude K > 22 
and the distribution of the Az/(1 + z) for these sources has 
values for the mean and a consistent with the ones derived 
for brighter sources (see right panel of Figure 1). The statis- 
tics at those faint magnitudes is not large enough to allow 
us to completely exclude systematics at z > 2 due to low 
signal-to-noise in the K-b&nd and specially in the Spitzer 
IRAC bands. However, it is worth noting that for this work 
we focus at z < 2, a range in which our photometric redshifts 
are more reliable, as shown in Figure 1. 

Finally, it is also worth noting the good accuracy of our 
photometric redshift estimates at high redshift, z > 4. In 
particular, in the field there is a faint galaxy which has been 
detected in the if-band {K A b= 23.17) with WFCAM and 
which has a very high spectroscopically confirmed redshift 
of z = 5.6 (Ouchi et al. 2005). For this galaxy our analysis 
yields a photometric redshift of 2 = 5.5±0.2, providing addi- 
tional confidence in the robustness of our redshift estimation 
procedure out to very high redshifts. 

4.2 The K-z relation 

Figure 2 shows the K — z relation for the sources in the 
final UDS-galaxy sample. This plot illustrates the coverage 
of the K — z plane provided by the sample, and also pro- 
vides a useful sanity check that not only the overall redshift 
distribution, but also the photometric redshift values for in- 
dividual sources are reasonable. 

The dashed line represents the empirical K — z relation 
for radio galaxies as obtained by Willott et al. (2003) , which 
approximately coincides with the passive evolution of a 2>L* 
starburst formed at redshift Zf — 10. The solid line cor- 
responds to the passive evolution of a 10L* starburst with 
the same formation redshift. To help minimise catastrophic 
errors in the photometric redshift estimation procedure, we 
introduced a prior which required that the redshift of the 
sources should not exceed the K — z relation of a passively 
evolving 10L* (i.e. the solid line). In fact this prior was only 
required for 12 galaxies, and it can be seen that the distribu- 
tion of galaxies on the K — z plane appears reasonable, with 
no significant excess of galaxies towards the bright limit. 

4.3 The BzK diagram 

Another consistency check on our photometric redshifts is 
offered by the BzK diagram. As expected, all the sources 



with photometric z > 1.4 lie in the regions with BzK = 
(z'-K)-(B-z') > -0.2 and BzK < -0.2 n (z'-K) > 2.5 
discussed by Daddi et al. (2004). In fact due to sheer number 
of sources in the UDS EDR sample, the BzK diagram dis- 
plays a number of interesting features of relevance to galaxy 
evolution, some of which have not been evident in any pre- 
vious studies. The BzK diagram, along with a discussion of 
these features, is discussed in detail by Lane et al. (2007) . 



4.4 Comparison with previous studies 

The redshift distribution of our new UDS galaxy sample is 
shown in Figure 3. The left-hand panel shows (as the solid 
histogram) the redshift distribution for the subset of sources 
with Kab< 21.9 (which corresponds to Kve ga < 20) to fa- 
cilitate comparison with the distribution obtained for the 
complete sample of 545 spectroscopically-confirmed sources 
(down to Kvega= 20) in the K20 survey (Mignoli et al. 
2005). The right-hand panel instead shows the redshift dis- 
tribution of the UDS galaxy sample down to Kab< 22.5 this 
time compared with the equivalent result obtained in the 
GOODS Chandra Deep Field South (CDFS) by Caputi et al. 
(2006) down to the same magnitude limit. Nearly half of the 
sources in the GOODS/CDFS have spectroscopic redshifts 
and the remainder have photometric redshifts computed by 
using 12 broad-band photometry from 0.4pim to 8/im (Cira- 
suolo et al. - in preparation). Both the distributions from 
K20 survey and GOODS/CDFS have been re-normalised to 
the same area as the UDS-galaxy sample. 

We find a very good agreement between the redshift dis- 
tribution we obtain for sources in the UDS galaxy sample 
and the ones displayed by the smaller samples which, how- 
ever, do obviously possess a much higher fraction of spectro- 
scopically confirmed redshifts. It is also evident from Figure 
3 that the surface density of galaxies in the present sam- 
ple is in agreement with those derived from these previous 
studies. The consistency of the redshift distributions and of 
the surface densities at different magnitude limits provides 
further confidence in the reliability of our sample and our 
photometric redshift estimates. 



5 LUMINOSITY FUNCTION 

In this Section we exploit the large area and wide redshift 
coverage of the UDS galaxy sample to derive the evolution 
of the near-infrared LF from the local Universe up to z > 2. 
The rest-frame K-band magnitudes for sources in the sam- 
ple have been computed by using the closest observed band 
depending on the redshift of the source. This method in fact 
minimises the uncertainties related to the fe-corrections, even 
though in the K-band these corrections have a small vari- 
ance compared to the optical wavelengths (Cowie et al 1994; 
Poggianti 1997). In our case the availability of the IRAC 
3.6/im and 4.5/im data enables us to directly trace the rest- 
frame K-band up to z ~ 1 for more than half of the sam- 
ple. However, as a test we also compared the absolute Mk 
magnitude obtained by using IRAC data with that obtained 
by directly extrapolating the observed K-band magnitude. 
The distribution of the difference of the magnitudes com- 
puted with these two different methods is centred around 
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Figure 3. The estimated redshift distribution of the new UDS galaxy sample, compared with previously published redshift distributions 
from near-infrared samples of comparable depth but much smaller area. The solid histogram shows the distribution of photometric 
redshifts for sources in the UDS galaxy sample with K^s< 21.9 and Kab< 22.5 in the left and right panels respectively. For comparison, 
the redshift distributions of galaxies from the K20 survey (Cimatti et al. 2002; Mignoli et al. 2005) and GOODS/CDFS (Caputi ct al. 
2006; Cirasuolo et al. in preparation) are also plotted as shaded histograms in the left and right panels respectively, after normalisation 
to the same area as the UDS galaxy sample. 



zero with rms a ~ 0.35. Therefore, these two methods pro- 
vide consistent results and for our dataset we used the Mis- 
estimated by using the IRAC bands when available. 

Figure 4 shows the rest-frame A-band luminosity func- 
tion (LF) in the redshift range 0.25 < z < 2.25 computed in 
six redshift bins by using both the classical 1/V ma x method 
(Schmidt 1968) and the maximum likelihood analysis (Mar- 
shall et al 1983). The 1/V ma x method does not require any 
assumption regarding the shape of the LF or parameter de- 
pendence but can be biased by strong density inhomogene- 
ity in the field. The solid dots in Figure 4 represent the 
LF obtained with this method. Only the bins above the 
luminosity-completeness limit at each redshift bin are shown 
and the error bars on each point correspond to the Poisso- 
nian errors. 

The solid line in Figure 4 shows the LF obtained in each 
redshift bin with the maximum likelihood method. This is a 
parametric method, for which we chose a Schechter function 
to describe the shape of the LF: 



Table 1. Best fit parameters for the Schechter LF. 



4>{M) = 0.41n(10)</>*10~ 04AM(Q+1) exp(-10~ 



(1) 



with AM = M K - M* K . The overall normalisation (<£*) 
has been obtained by matching the number of observed 
galaxies in each redshift bin. The best fit values for the pa- 
rameters a, Mk and 0* in each redshift bin are shown in 
Table 1. 

In the determination of the LF both in the 1/V ma x 
method and maximum likelihood analysis a correction factor 
has been applied to take into account the incompleteness ef- 
fects. In fact, as described in Foucaud et al. (2007), the UDS 
EDR sample is nearly 100% complete at Aab< 21.9 and the 
completeness drops to ~ 70% at Aab = 22.5. 



z 
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M K 


0*(1O- 3 Mpc~ 3 ) 


0.25 - 


0.75 


-0.99 ±0.04 


-22.86 ± 0.08 


2.9 ±0.2 


0.75 - 


1.00 


-1.00 ±0.08 


-22.86 ± 0.11 


3.4 ±0.3 


1.00 - 


1.25 


-0.94 ±0.13 


-22.93 ± 0.14 


2.8 ±0.4 


1.25 - 


1.50 


-0.92 ±0.18 


-23.03 ± 0.16 


2.6 ±0.3 


1.50 - 


1.75 


-1.00 fixed 


-23.23 ± 0.09 


1.7 ±0.2 


1.75 - 


2.25 


-1.00 fixed 


-23.58 ±0.13 


0.6 ±0.3 



Due to the limiting magnitude A"ab< 22.5 our ability 
to trace the faint end of the LF decreases with redshift. How- 
ever, up to z ~ 1.5 the likelihood analysis suggests that the 
value of the faint end slope a does not change with redshift 
and stays consistent with the local value of a ~ 1. The actual 
data do not allow us any estimate of the slope for z > 1.5. 
Therefore, in the likelihood analysis at redshifts larger than 
z = 1.5 we decided to fix the value of the faint end slope to 
the mean value obtained at lower redshifts (a — 1). This is 
consistent with the results of deeper surveys (e.g. Caputi et 
al. 2006). Recently, by exploiting a sample of ~ 300 objects 
selected from the Hubble deep field south (HDF-S) down to 
Aab< 24.9 Saracco et al. (2006) found the value of the slope 
of the LF to be constant with redshift and consistent with 
the local value (a ~ 1) up to z ~ 3. 

The agreement between the estimation of the LF ob- 
tained with the two independent methods (1/Vmax and max- 
imum likelihood) is evident in Figure 4. We also found our 
LF to be consistent with the results derived by Pozzetti et 
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Figure 4. Rest-frame i\-band luminosity function in six redshift bins. The solid dots show the LF obtained from the 1/Vmax method 
for sources in the UDS sample with Kp^B< 22.5. The open symbols correspond to the LF obtained from the K20 survey (Pozzetti et 
al. 2003) in three redshift bins ([0.2-0.65], [0.75-1.3], [1.3-1.9]), while the solid line is the best fit Schechter function obtained from the 
likelihood analysis. The LF obtained from the GOODS/CDFS sample (Caputi et al. 2006) is indicated by the long dashed line. Also 



shown, for comparison, are the LFs at z 
respectively. 



: (dashed line) and z = 1 (dotted line) from Kochanek et al. (2001) and Drory et al. (2003), 



al. (2003) from the K20 survey and Caputi et al. (2006) 
from GOODS/CDFS. It is worth noting that in this work 
we have nearly a factor of 10 more objects compared to the 
GOODS/CDFS sample and a factor ~ 40 more than the K20 
survey. Due to the large area we also have enough statistics 



to be able to properly trace the bright end of the LF over 
our full redshift range with negligible cosmic variance. 

The comparison with the local LF obtained by 
Kochanek et al (2001) from the Two Micron All Sky Survey 
(2MASS) shown in Figure 4 confirms the substantial bright- 
ening of the characteristic magnitude M% with redshift. In 
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agreement with previous studies (Pozzetti et al 2003; Drory 
et al. 2003; Feulner et al. 2003; Caputi et al. 2006; Saracco 
et al. 2006) we find M£ to be brighter by ~ 1 magnitude at 
z ~ 2 compared to the local Universe. The lack of data at 
the faint end of the LF does not allow us to precisely trace 
the density evolution of the LF. However, up to z ~ 1.75 
where the likelihood estimate of the best fit parameters for 
the LF is still reliable (see Table 1) we observe a decrease in 
the space density by a factor ~ 2 compared to z = 0. This 
decrement increases up to a factor ~ 4.5 if the highest red- 
shift bin is considered. This is in good agreement with the re- 
sults obtained by Caputi et al. (2006) in the GOODS/CDFS 
as also shown in Figure 4. However, the much large statis- 
tics obtained with the current work allows for substantial 
improvement in determining the LF, especially at its bright 
end. A more detailed likelihood analysis of the LF, including 
the parameterisation of the evolution and the determination 
of the faint end slope up to the highest redshift bin will be 
possible with the next data release of the UDS which will 
reach Kab^ 23.5. 



6 THE COLOUR BIMODALITY 

Following the above analysis of the evolution with redshift 
of the population of K-b&nd selected galaxies it is also in- 
teresting to investigate the role of the different galaxy pop- 
ulations and how they contribute to the global evolution of 
the JT-band LF. For this purpose we analysed the rest-frame 
(U — B) colour of sources in the UDS-galaxy sample; accord- 
ing to the well-studied local colour-magnitude relation (e.g. 
Visvanathan & Sandage 1977; Bower et al. 1992; Strateva 
et al. 2001) the red and blue galaxies show a bimodal dis- 
tribution which allows one to easily separate the population 
of red/early type galaxies from the blue/late type systems. 

Figure 5 shows the rest-frame (U — B) colour versus 
absolute magnitude in the B-band in the same six redshift 
intervals used for the LF. The rest- frame U and B band mag- 
nitudes have been computed by using the best fit template, 
according to the redshift of the source. The small insets in 
each panel show the distribution of the rest-frame (U — B) 
colour in each redshift bin. In agreement with studies in the 
local Universe (Strateva et al. 2001; Hogg et al. 2002; Blan- 
ton et al. 2003a,b; Baldry et al. 2004) and at z ~ 1 (Bell et 
al. 2004a; Weiner et al. 2005; Willmer et al. 2006; Franzetti 
et al. 2007) we find that the colour-magnitude relation ex- 
hibits a clear bimodality up to z ~ 1. 

The data also show the luminosity-dependent behaviour 
of the bimodality (Aragon-Salamanca et al. 1993; Stanford 
et al. 1995, 1998; Ellis et al. 1997; Kodama et al. 1998; 
van Dokkum et al. 2000) with the brighter sources display- 
ing redder colours. In Figure 5 the diagonal line represents 
the colour-magnitude relation derived by van Dokkum et al. 
(2000) for red galaxies in distant clusters adding an offset of 
-0.25 magnitudes to separate the red and blue populations 
(see Willmer et al. 2006): 

(J7-S) = -0.032(M B + 21.52) + 1.4 (2) 

converted to the AB magnitude system. As shown in Figure 
5, this luminosity-dependent colour cut is in good agreement 
with the observed bimodality present in the UDS galaxy 
sample, at least up to z ~ 1.5. 



However, it is interesting to note that the strength of the 
colour bimodality fades with increasing redshift and nearly 
disappears at z > 1.5. The extend to which the apparent 
disappearance of the colour bimodality is due to photomet- 
ric errors will be revealed as the UDS data progressively 
improve. A more detailed study on the evolution of the bi- 
modality with redshift is beyond the scope of the present 
paper and will be addressed in detail by Eales et al. (2007 - 
in preparation). 

However, it is worth noting the difficulty of the broad 
band rest-frame colours to distinguish between the different 
populations of galaxies that show the same optical and near- 
infrared colours, i.e. the old passive evolving sources and the 
dusty star-forming one. As well demonstrated in a recent 
paper by Stern et al. (2006) dusty starbursts and passively 
evolving galaxies are virtually indistinguishable in the broad 
band photometry at wavelengths < 10/im. 

In the local Universe, the population of red galaxies 
is mostly constituted by morphologically early type objects 
(Sandage & Visvanathan 1978; Bower et al. 1992; Terlevich 
et al. 2001) and the mixing between the populations of early 
and late type objects is relatively low. Recent studies on the 
SDSS survey by Strateva et al. (2001) found only ~ 20% of 
late type galaxies to have red colours compatible with the 
red sequence. 

This low contamination by the late type galaxies in the 
red sequence has also been observed at intermediate redshift 
( z < 1) both in the field (Kodama et al. 1999; Bell et al. 
2004b) and in clusters (Couch et al. 1998; van Dokkum et 
al. 2000; van Dokkum & Franx 2001). These results are in 
good agreement with the well defined bimodality we observe 
at z < 1 in the UDS-galaxy sample (see Figure 5). 

However, the contamination and the mixing between 
early and late type galaxies in the colour-magnitude dia- 
gram increases with redshift as clearly shown in Figure 5. 
Recently, Franzetti et al. (2007) exploiting data at z ~ 1 
from the VIMOS-VLT Deep Survey (VVDS) found that even 
though the red peak of the colour distribution is mainly pop- 
ulated by early type galaxies and the blue peak by late type 
ones, there is a strong contamination between the two pop- 
ulations. In particular, they found that in the population of 
red galaxies only ~ 65% are early type galaxies, while the 
remaining ~ 35% are late type objects. 



7 THE LUMINOSITY FUNCTION 
SUBDIVIDED BY COLOUR TYPE 

As discussed in the previous Section and in Eales et al. 
(2007 - in prep.), classification by colour into red and blue 
sources provides an efficient means of identifying the popu- 
lation of early and late type objects at z < 1. However, as 
also outlined by Eales et al. the difficulties in successfully 
exploiting this colour selection increase with increasing red- 
shift as demonstrated by the progressively disappearance of 
the colour bimodality at z > 1. Therefore, in the following 
Section we will study the evolution with redshift of the pop- 
ulations of red and blue galaxies bearing in mind that this 
roughly corresponds to the evolution of the early and late 
type objects only at low and intermediate redshift. How- 
ever, it is worth remembering that even at these redshifts 
the cross-contamination of each population can be relatively 
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Figure 5. Rest-frame (U — B) colour versus absolute magnitude in the B band in six redshift intervals for sources in the UDS galaxy 
sample. The small insets show the distribution of the rest-frame (U — B) colour in each redshift bin. The diagonal line represents the 
luminosity-dependent colour cut (see Equation 2) adopted by van Dokkum et al. (2000) to divide the populations of red and blue galaxies. 



high (~ 20 - 30% ; Strateva et al. 2001; Ilbert et al. 2006; 
Franzetti et al. 2007). 

As a first step we sub-divided the UDS galaxy sam- 
ple adopting the luminosity-dependent colour cut shown by 
Equation 2 (van Dokkum et al. 2000). The resulting LFs 
for the red and blue galaxies are shown in Figure 6. As for 
the LF of the total sample of A"-band selected galaxies (see 
Section 5) we used both the 1/V ma x and maximum likeli- 
hood methods. The best fit parameters for the likelihood are 



shown in Table 2. Due to the limiting magnitude Kab< 22.5 
the determination of the best fit parameters in the redshift 
bin 1.75 < z < 2.25 is particularly challenging. Therefore, 
in the following discussion we consider only the redshift bins 
in the range 0.25 < z < 1.75 and simply comment on the 
behaviour in the highest z bin. 

The limiting magnitude Kab< 22.5 also does not allow 
us to obtain a precise determination of the faint end of the 
LF, as already discussed in Section 5. Therefore, the maxi- 
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Figure 6. Rest-frame l\-band LF for sources in the UDS-galaxy sample sub-divided by colour according to Equation 2. The solid dots 
represent the LF of the population of red galaxies while the open triangles show the LF of the blue objects. The dotted and long dashed 
lines show the LF obtained in each redshift bin by using the maximum likelihood analysis for the red and blue population respectively 
(see Table 2). The solid line is the LF of the whole i-f-band selected sample in each redshift bin as plotted in Figure 4 and tabulated in 
Table 1. For reference we also indicate by a dashed line the local if -band LF from Kochanek et al. (2001). 



mum likelihood analysis has been performed both leaving a 
as a free parameter as well as fixing this value to the one ob- 
tained in the lower redshift bins, where our data better trace 
the faint-end slope. The mean value for the parameter a in 
the redshift range 0.25 < z < 1.5 is -1.35 and -0.1 for the 
blue and red objects respectively. Therefore, at z > 1.5 we 
fixed the slope of the LF at these values. As is clearly seen 



in Figure 6, the population of blue galaxies exhibits a much 
steeper slope compared to the red galaxies up to z ~ 1.5 in 
agreement with previous works (Lilly et al. 1995; Blanton 
et al. 2001; Pozzetti et al. 2003; Giallongo et al. 2005; Ilbert 
2006). 

The redshift evolution of the blue and red populations 
also appears to differ. Even though both the red and blue 



11 



Table 2. Best fit parameters for the Schechter LF of red and blue 
galaxies. 



Ml 



</>*(icr 3 Mpc" 3 ) 



0.25 - 0.75 
0.75 - 1.00 
1.00 - 1.25 
1.25 - 1.50 
1.50 - 1.75 
1.75 - 2.25 



Blue objects 
-1.19 ±0.05 -22.71 ±0.11 



-1.44 ±0.15 
-1.38 ±0.19 
-1.43 ±0.24 

— 1.35 fixed 

— 1.35 fixed 



-23.05 ±0.18 
-23.17 ±0.21 
-23.40 ±0.25 
-23.45 ±0.15 
-23.75 ±0.20 



1.5 ±0.2 
1.3 ±0.2 
1.3 ±0.3 
1.0 ±0.4 
0.9 ±0.3 
0.25 ±0.3 



Red objects 



0.25 - 0.75 
0.75 - 1.00 
1.00 - 1.25 
1.25 - 1.50 
1.50 - 1.75 
1.75 - 2.25 



-0.12 ±0.1 
-0.09 ±0.2 
-0.05 ±0.4 
-0.14 ±0.5 
-0.1 fixed 
-0.1 fixed 



-22.36 ±0.07 
-23.48 ±0.13 
-22.65 ±0.17 
-22.65 ±0.24 
-22.81 ±0.10 
-23.04 ±0.21 



2.3 ±0.2 
2.0 ±0.3 
1.0 ±0.3 
1.0 ±0.4 
0.7 ±0.3 
0.2 ±0.4 
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Figure 7. Space density as a function of redshift for bright 
sources with Mk < —23 classified as red (solid dots) and blue 
(open triangles) according to their rest-frame colours (see Equa- 
tion 2). The black solid squares represent the space density of all 
sources with Mk < —23. 



sources experience a brightening of the characteristic magni- 
tude M K , the blue objects show a slightly stronger evolution 
with AM K ~ -0.7 over the redshift range 0.25 < z < 1.75 
compared with AA1 K ~ —0.5 of the red population. Accord- 
ing to the results of the maximum likelihood analysis, both 
populations show a small density evolution at z < 1 while at 
higher redshifts the decrease in number density is stronger 
for red objects compared with the blue ones. However, it is 
worth remembering that the uncertainties in the classifica- 
tion of the two populations based purely on the rest-frame 
colour cut expressed by Equation 2 can introduce severe 



bias, specially at higher redshift where the bimodality fades 
away. 

The inter-connection between the luminosity and den- 
sity evolution coupled with the different faint-end slope of 
the LFs for the two populations result in the predominance 
of the blue objects at faint magnitudes (Mk <; —22) over 
the whole redshift range explored in this work. On the other 
hand, the bright end of the LF is dominated by red objects 
(see Figure 6) but only up to z ~ 1. In fact, as shown in Fig- 
ure 7 the space density of bright sources with Mk < —23 is 
nearly constant for red objects in the range 0.25 < z < 1.5, 
while it increases by a factor ~ 2 for the blue ones over the 
same redshift interval. The result is that the population of 
the brightest objects mostly consists of red sources at low 
and intermediate redshift, while at z > 1 there is an inver- 
sion in the proportion with a slight predominance of blue 
galaxies. 

It is also worth noticing that the mild density evolution 
of the red galaxies between z = and z ~ 1 is in good agree- 
ment with the recent results by Yamada et al. (2005). These 
authors found that the number density of old bright passive 
galaxies is decreased by only a small factor (less than 2) 
compared with the local Universe. This is inconsistent with 
the results of Bell et al. (2004a) which instead suggested a 
decrease in the population of red galaxies by nearly a fac- 
tor 10 over the same redshift range. However, as shown by 
Bundy et al. (2006) this decrement can be due to a selection 
effect. In fact, these authors pointed out how the J?-band 
limit - the same band used by Bell et al. to select their 
sample - introduces a bias against red galaxies, especially 
at z > 1, and which introduces an artificially large decrease 
with redshift in the red population. 



8 DISCUSSION AND CONCLUSIONS 

We have presented new results on the evolving near-infrared 
galaxy luminosity function, based on the analysis of a 
large sample of 7f-band selected galaxies drawn from the 
Early Data Release of the UKIDSS Ultra Deep Survey. 
The complete sample includes ~ 22, 000 sources down to 
Kab< 22.5 over an area of 0.6 square degrees. Exploiting 
the complementary optical and mid-IR data provided by 
the Subaru/XMM-Newton Deep Survey and Spitzer-SWIKE 
Survey we have obtained reliable photometric redshift es- 
timates for the sources in the sample. This allowed us to 
study the evolution of the Jf-band LF in the redshift range 
0.25 < z < 2.25. In agreement with previous studies, we 
find that the characteristic luminosity of the near-infrared 
luminosity function exhibits a substantial brightening of ~ 1 
magnitude between the local Universe and z ~ 2, while the 
total density decreases by a factor ~ 2 over the same redshift 
interval. 

We have also used the rest frame (U — B) colour to 
separate the sample into red and blue galaxies. We find the 
classical bimodal colour distribution to be clearly present 
at z < 1. However, the strength of the colour bimodality 
is found to decrease with increasing redshift and virtually 
disappears by z > 1.5. The progressively improvement of 
the UDS data will assess how the errors on the photometry 
affect this result. 

By using the colour - magnitude relation to split the 
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Figure 8. The solid dots represent the evolution of the char- 
acteristic luminosity with rcdshift derived from the likeli- 
hood analysis relative to the one in the local Universe (M^(z = 
0) = 22.26) as obtained by Kochanck ct al. (2001). The dotted 
line is the evolution of M* K derived by Caputi et al. (2006) from 
GOOD/CDFS. The solid and dashed lines represent the evolution 
with redshift of the _K"-band absolute magnitude derived by a sin- 
gle burst model from Bruzual & Chariot (2003) with formation 
rcdshifts z = 3 and z = 6, respectively. 



sample we have analysed the cosmological evolution of the 
red and blue galaxy populations. The blue objects clearly 
dominate the faint end of the LF over the whole redshift 
range explored in this work. On the other hand, up to z ~ 1 
the brightest sources are mainly provided by red objects. In 
fact, the space density of red sources with Mk < —23 is ap- 
proximately constant with redshift while the blue galaxies, 
with the same bright absolute magnitude, exhibit a strong 
increase in their space density with redshift. Therefore, this 
different evolutionary behaviour of the two populations de- 
termines the preponderance of blue objects in the bright-end 
of the luminosity function at redshifts z > 1. 

Compared with previous studies, the unique combina- 
tion of large area and depth of the UDS minimises the effect 
of the cosmic variance and has led to a substantial improve- 
ment in the determination of the LF. In particular, the large 
area coverage has allowed us to accurately trace the evolu- 
tion of the brightest sources at z < 2 with a high level of 
statistical significance. As shown in Figure 7, the space den- 
sity of the total population of bright sources {Mk < —23) 
is roughly constant over the redshift range 0.5 < z < 1.5. 
This strongly suggests the vast majority of the present-day 
bright elliptical galaxies were already in place at z ~ 1.5. 

On the other hand, our analysis has also shown a change 
in the fraction of red and blue galaxies with redshift. The 
majority of bright sources at z < 1 exhibit red colours con- 
sistent with having an old stellar population formed at high 
redshift. However, the increase of the fraction of blue galax- 
ies with increasing redshift - indicative of an increase of 



the ongoing star formation - suggests the major formation 
epoch of the these bright sources to be at z > 1.5. 

However, it is worth noting that the colour classifica- 
tion based on the rest-frame (U — B) colour is not very ac- 
curate for distinguishing between the different populations 
of galaxies that show the same optical colours, i.e. the old 
passive evolving sources and the dusty star-forming one. In 
fact, studies of the population of Extremely Red Objects 
(EROs) at redshifts 1 < z < 2 revealed this population to 
be composed by both old passive evolving galaxies and dusty 
star-bursts, roughly in the same proportion (Cimatti et al. 
2002; Smail et al. 2002; Yan & Thompson 2003; Cimatti et 
al. 2003; Gilbank et al., 2003; Moustakas et al. 2004; Simp- 
son et al. 2006b). 

On the other hand, not all the blue objects are experi- 
encing their major event of star formation. In fact, it only 
requires less than 10% of the total mass of an early type 
galaxy to be involved in a recent episode of star formation 
to completely alter the colour of the objects and show it as 
a blue galaxy (Zepf 1997). This is because even a relatively 
small number of bright young stars dominates the emission 
at optical wavelengths. Moreover, a recent study from the 
VVDS in the redshift range 0.4 < z < 0.8 by Ilbert et al 
(2006) found that ~ 30% of bulge dominated sources exhibit 
blue colours, indistinguishable from the late type population 
(see also Im et al. 2001; Menanteau et al. 2004; Cross et al. 
2004). 

It is also interesting to notice that the brightening of 
the characteristic luminosity (M K ) of the whole Tf-band 
selected LF is consistent with the passive evolution of a 
single stellar population. Figure 8 shows the evolution of 
M K derived from our likelihood analysis relative to M K at 
z=0 (Kochanek et al. 2001). Within the errors, our data 
agree remarkably well with the evolution of the Tf-band ab- 
solute magnitude displayed by a single burst model from 
Bruzual & Chariot (2003) with formation redshifts in the 
range 3 < z < 6. This result, coupled with the constant 
space density of bright sources up to z ~ 1.5, suggests the 
major assembly epoch of bright sources to be at z > 2, 
even though spread over a large range in redshift. In fact, 
as shown from Figure 7, nearly 40% of the bright sources at 
z ~ 1.5 exhibit red colours. Assuming that half of them are 
dusty starbursts - as suggested by studies of ERO popula- 
tions (e.g. Cimatti et al. 2002) - we end up inferring that 
at least 20% of sources are old passively evolving galaxies 
at these redshift. Therefore, this implies that a reasonable 
fraction of bright sources at z ~ 1.5 was assembled at very 
high redshift with very efficient star formation. 
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